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We present a study of the decays B°' + — > J/ipujK ' + using 383 x 10 6 BB events obtained with the 
BABAR detector at PEP-II. We observe y(3940) -> J/ipu), with mass 3914.6±3;4(s£ai) ± 2.0(syst) 
MeV/c 2 , and width 34±g 2 (siai) ± 5(syst) MeV. The ratio of B° and B + decay to YK is 
0.27lQ'23(stai)^g'gj(s?/st), and the relevant B° and B + branching fractions are reported. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The BELLE Collaboration has reported evidence for 
the A(3940) 0, the F(3940) 0, and the Z(3930) [3. 
The mass and width values are the same within error, 
the states have positive C parity, and spin-parity (J p ) 
2 + is favored for the Z, which may then be the first 
radial excitation of the % C 2(3556), i.e., a charmonium 
state. The mass and width consistency with the X and 

Y suggests the possibility that these may be the Z in 
different production contexts. The Z was found in two- 
photon production of DD, so that it may be a charmo- 
nium state. The X was observed in e + e~ — ► J/ipX, and 
decays mainly to D*D, suggesting a charmonium inter- 
pretation. In contrast, the Y was found in B — > YK, 

Y — > J/ipioj, which is OZI suppressed for a charmo- 
nium state [4J. Also, an analysis of B — > KDD and 
B — > KD*D [5| shows no evidence for the Y (nor for the 
X or Z), although ^(3770) -> £>S and X(3872) -> £>*S 
are observed. Other possibilities for the nature of this 
state, already suggested for the X(3872), include a hy- 
brid charmonium-gluon bound state,cco In Ml, a molec- 
ular state of a cc(uu + dd) system [8j, |9j, [lfj, |ll|, [l2| , or 
a multiquark state [13J. The 5— wave molecule model [9J 
predicts a very small B°/B+ ratio for B -> KX(3872). 
The previous low value has been confirmed [14[ , although 
the uncertainties are still large, so that a measurement 
of this ratio for the F(3940) may be important to an 
understanding of this state. 

In this Letter, we examine the decays B 0,+ — ► 
J/iPir+Tr-ir°K°' + [15], with 71 



7r mass in 



the 



gion. We confirm the F(3940), improve the precision 
of the mass and width significantly, and measure the 
(B°/B + ) production ratio for the first time. Branch- 
ing fraction values for B — > YK, Y — > J/ipu>, and for 
B — ► J/ipujK are obtained for 2?° and B + decay sepa- 
rately; each is a first measurement. 

The data were collected with the BaBar detector [16( 
at the PEP-II asymmetric-energy e + e~ storage rings op- 
erating at the T(45) resonance. The integrated lumi- 
nosity for this analysis is 348 fb _1 . The decays B° :+ —* 
J/ipir + Tr~ir°K°' + are reconstructed as follows (TableHJ. 
A candidate J/ip — ► e + e~ (/i + /i~) decay has invariant 
mass in the J ftp mass region, and is then constrained 
to the nominal mass [rjj. A Kg candidate has 7r + 7r~ 
invariant mass in the Kg region. The J/tp and Kg dis- 
tributions from the B signal region show no significant 
background. A ir° candidate consists of a photon pair 
with invariant mass in the ir° region. After a ir° mass 
constraint, an uj —* 7r + 7r _ 7r° candidate has invariant mass 
in the u) region. We form a B + (B°) candidate by com- 



E* B 



bining J/ip, u and K + [18( (Kg) candidates. 

We define the B signal region using the center 
of mass (cm.) energy difference AE 
y/s/2, and the beam-energy substituted mass tties 
yf((s/2+pi ■ p B ) I Ei) 2 - p£ [li], where [E^pi) is the 
initial state four-momentum vector in the laboratory 
frame (l.f.); y/s is the cm. energy, E B is the B me- 
son energy in the cm., and pb is its l.f. momentum. 
Signal events have AE ~ zero and tties ~ tub, 12% of 
the events have multiple candidates, and for these the 
combination with the smallest \AE\ is chosen. 

The selection criteria were established by optimizing 
signal-to-background ratio using Monte Carlo (MC) sim- 
ulated signal events, B — > YK, Y — > J/ipw, and back- 
ground BB and e + e~ — ► qq (q = u, d, s, c) events. 

TABLE I: Principal criteria used to select B candidates. 



Selection Category 


Criterion 


J/tp —* M + M _ mass (GeV/c 2 ) 


3.06 < m MM < 3.14 


J/tp — > e + e _ mass (GeV/c 2 ) 


2.95 < m ee < 3.14 


Ks mass (GeV/c 2 ) 


0.472 < mnn < 0.522 


7T° mass (GeV/c 2 ) 


0.115 < m 77 < 0.150 


uj signal region (GeV/c 2 ) (B + ) 


0.7695 < m 3 n < 0.7965 


uj signal region (GeV/c 2 ) (B°) 


0.7605 < m 3 n < 0.8055 


AE (GeV) (B+) 


\AE\ < 0.020 


AE (GeV) (B°) 


\AE\ < 0.015 


m ES (GeV/c 2 ) 


5.274 < tubs < 5.284 


B helicity angle 9b 


|cos6» s | < 0.9 


Photon helicity angle Q 1 


cos# 7 < 0.95 


tp (2S) veto ( GeV/c 2 ) 


3.661 < Mj/^ v < 3.711 



The cos6*b distribution (8b is the cm. polar angle of 
the B) is proportional to sin 9b', since e + e~ — ► qq events 
peak toward ±1, we require | cos 9b | < 0.9. The variable 
cos 9j is the normalized dot product between the higher 
momentum photon in the ir° rest frame (r.f.) and the l.f. 
direction of the ir°. For ir° decay this distribution is flat; 
background peaks at 1, hence we require cos# 7 < 0.95. 
Events from B -> ip (25) Kit , ip (2S) -> ir+ir-J/iP, are 
removed by the ip (25) veto. 

The 37T mass, rriESi an d AI? distributions are shown 
in Fig. [TJ where we apply all Table U criteria except the 
requirement on the variable plotted. We fit the Sir mass 
distributions with an w-meson Brcit-Wigncr (BW) line 
shape (nominal u) mass and width [l7|) convolved with 
a MC-determined triple-Gaussian resolution function as 
signal, and a quadratic background function. We fit the 
nigg distributions with a signal Gaussian with mass and 



width fixed from MC, and an ARGUS background func- 
tion [l9(, and fit the AE distributions with a double- 
Gaussian signal function determined from MC, and a lin- 
ear background function. 
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FIG. 1: (a)-(c) ((d)-(f)) The 3vr mass, m E s, and AE dis- 
tributions for the B + (B°) mode; solid (open) dots are for 
unweighted (weighted) events. The solid (dashed) curves rep- 
resent signal plus background (background). 

There is a large u> signal for the B + mode, and a 
smaller signal for B°; the tues and AE distributions 
exhibit clear B signals. We establish the correlation be- 
tween the ui and B signals with a projection procedure 
based on the to decay angular distribution. The helicity 
angle, Oh, is the angle between the 7r + and 7r° directions 
in the 7r + n~ r.f.. The cos 6^ distribution is proportional 
to sin 2 Oh, and the ui signal is projected by giving the i th 
event weight Wi = |(1 — 3cos 2 #^). The effect is shown 
in FigQ] For the B + mode, the omega signal survives, 
and background is removed. For the B° mode the ef- 
fect is qualitatively similar. Confirmation is obtained 
from a fit to the 3-7T mass distribution in each interval. 
We conclude that there is one-to-one correspondence be- 
tween the uj and .B-meson signals in tues & n d AE, and 
that, at the present level of statistics, the 3ir system 
in the uj mass region results entirely from m decay for 
B —> J/-i/j7r + ir~ir K. The uj — ities (or AE) signal corre- 
lation is important to an analysis of the J/ipuj threshold 
mass region. Near threshold, the 3tt mass distribution 
above the u> mass is limited in range and distorted in 
shape. The ities distribution is not affected, and so we 
use rriES fits to extract the J/ipu) mass distribution. 

For each B decay mode, the ities distribution in each 
interval of J / ipir + ir~ ir° invariant mass is fitted to ex- 
tract the J/4>UJ signal. The tues signal, and ARGUS 



background, functions are those of Fig. [TJ the fits use a 
binned Poisson likelihood function with signal and back- 
ground normalizations free [201 ] . All fits converge prop- 
erly and provide good descriptions of the data. From 
threshold to 4 GeV/c 2 , the J/tpuj mass resolution varies 
from 5 — 8 MeV/c 2 , and so in this region the spectrum 
is investigated in 11 intervals of 10MeV/c 2 starting at 
3.8725 GeV/c 2 . At higher mass, there is no evidence of 
narrow structure, and we show the results in 50 MeV/c 2 
intervals. In Fig. [2j there is a clear enhancement near 
threshold for B + decay, while at higher mass no struc- 
ture is apparent. The total B + (B°) signal in Fig. [2] is 
236±^ (32±|) events of which 109±i| (Wt 7 6 ) have J/ijjuj 
mass less than 4 GeV/c 2 (statistical errors only). 

We correct the mass distributions of Fig. [2] for effi- 
ciency and resolution. In the MC simulation of the Y 
signal, we assume phase space decays of B — ► YK and 
Y — ► J/ij)U), but use the correct angular distribution for 
uj decay. Initially we used a relativistic S-w&ve BW line 
shape with M(Y) = 3.940 GeV/c 2 and T(Y) = 0.06 
GeV [2|. Mass resolution effects result in a net flow 
of events away from the peak mass value. For a given 
mass interval we define acceptance as the ratio of events 
reconstructed in that interval to events generated in the 
interval; this accounts for efficiency and resolution effects. 
The acceptance-corrected spectrum is fit to a relativistic 
BW line shape without convolving resolution, since the 
acceptance correction takes this into account. We obtain 
values of M(Y) and T(Y) which are smaller than in the 
initial simulation, and so generate new MC samples with 
the new values in order to correctly reproduce resolu- 
tion effects. This iterative procedure converges quickly, 
and the acceptance results in Figs.[2]Jc), (d) are obtained 
with M(Y) = 3.915 GeV/c 2 and T(Y) = 0.02 GeV. The 
dip at ~ 3.91 GeV/c 2 is due to net flow of events away 
from the resonance maximum because of mass resolution. 
At lower mass, the acceptance is slightly lower than at 
higher mass because of the proximity to threshold. Al- 
though the acceptance variation in the Y signal region 
is significant, the effect on the Y fit parameters, and on 
the corrected number of signal events, is small because 
of the large statistical uncertainties on the data. 

The decrease in acceptance at high mass in Fig. [2fd) 
results from decreasing Kg l.f. momentum. The decay 
pion reconstruction probability decreases because its l.f. 
momentum is too small, or because the decay opening 
angle is so large that the pion does not intersect enough 
detector planes. Fig.[3]shows the corrected mass distribu- 
tions. Below ~ 4 GeV/c 2 we correct interval-by-interval, 
while for higher mass we use a linear fit to the J/iJjuj 
mass dependence. The B° data are corrected for K L and 
K° s -> 7T 7r° decays. 

We associate the near-threshold enhancement in 
Fig. (Ha) with Y production [2J, and obtain the mass, 
width and decay rate from y 2 fits. The fit function 
consists of a relativistic S— wave BW describing the Y 
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FIG. 2: The J/ipuj mass distribution from the Toss fits for 
(a) B + , and (b) B° decay. The acceptance as a function of 
J/ipuj mass (c) for the B + , and (d) for the B° mode. 



and a Gaussian nonresonant contribution. The cor- 
rected B + and B° distributions are fitted simultane- 
ously, with mass, width and Gaussian parameters as com- 
mon free parameters. The fit describes the data well 
(x 2 /NDF = 45/44, NDF=number of degrees of free- 
dom ), as shown in Fig. [3l In Fig. Ela), the acceptance- 
corrected number of events with J /ipuj mass less than 
3.98 GeV/c 2 is 2140 ± 290(siai), while for the Gaussian 
it is 420 ± 90(stat). Our average efficiency of ~ 5% im- 
plies that a background fluctuation of ~ 19 standard de- 
viations would be required to describe the near-threshold 
enhancement. This occurrence has negligible probability, 
and so we have instead a clear observation of the Y(3940). 
The simultaneous fit yields a Y signal of 1980t37g(stai) 
events (i.e. magnitude 5.2 standard deviations) for B + , 
and 527±Hl(stat) for B°. 

Since the acceptance-correction procedure may depend 
on the input MC Y(3940) line shape, we combine the 
first 11 mass intervals for data and MC and make an 
overall efficiency correction. The results differ by 1.9%, 
and we incorporate this as a systematic error associ- 
ated with the MC line shape. Other systematic errors 
are estimated by repeating the entire process, separately 
varying by ±lcr the signal peak and width, and the 
ARGUS parameter, for the itles fits. The largest sys- 
tematic uncertainty contributions to the B + branching 
fraction are 5 — 6% due to the uncertainties in the sec- 
ondary branching fractions, tracking efficiency, and par- 
ticle identification. For B°, the largest contribution is 
10% due to rriES mass variation; secondary branching 



fractions, particle identification, tracking and K~s recon- 
struction efficiency contribute also. For both modes, 
there are uncertainties associated with the number of BB 
events produced, and with MC sample size. The prod- 
uct branching fraction for B + — ► YK + , Y — > J/tpuj is 
(4.9tli(stat) ± 0.5(syst)) x 10~ 5 , and that for B° -> 
YK°, Y -> J/ipw is (1.3±H(stat) ± 0.2{syst)) x 10~ 5 , 
with upper limit (95% C.L.) 3.9 x 10~ 5 for the lat- 
ter. The corresponding branching fractions for B — v 
J/ipujK are (3.5 ± 0.2(stat) ± OA(syst)) x 10~ 4 , and 
(3.1 ± 0.6(stat) ± 0.3(syst)) x 10" 4 , respectively. 

We define Ry and Rnr as the ratios between the 
number of B° and B + events (after all corrections) 
for the Y signal and for the nonresonant contribution, 
respectively. Simultaneous fits to Figs. EJa),(b) yield 
the values R Y = 0.27±j^(stat)±^ 4 (sys£) and r nr = 
0.97tof 2 (stat)to l{systy, the upper limit (95% C.L.) on 
Ry is 0.75. Although the uncertainty is large, the cen- 
tral value of Ry is smaller than expected from isospin 
conservation. In comparison, R is 0.865 ± 0.044 for 
B -> J/ipK [II and 0.81 ± 0.05(siai) ± O.Ol(syst) for 
B -* i){2S)K |2|. 
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FIG. 3: The corrected J/ipu) mass distribution for (a) B + and 
(b) B° decay. Each solid (dashed) curve represents the total 
fit function (the nonresonant contribution). 



The Y mass and width measurements are subject to 
additional systematic effects. When MC-generated signal 
events arc fitted using the input line shape with mass and 
width as free parameters, the fitted value of the mass is 
1.6 MeV/c 2 lower than the input value of 3.915 GeV/c 2 . 
This results from the limited 37r phase space near J /tpui 
threshold, and so we increase the fitted Y mass value by 
1.6 MeV/c 2 , and assign this as a systematic uncertainty. 



Also, we have used an S'-wave BW line shape to describe 
the Y. We repeat the fit using a P-wave line shape. The 
fitted mass value decreases by lMeV/c 2 , and the width 
increases by 5 MeV. We assign these as systematic uncer- 
tainties due to the choice of orbital angular momentum. 
Finally, a fit to the uncorrected distributions (Fig. \2§ 
yields a mass value f.4 MeV/c 2 larger, and a width 4 
MeV larger, than obtained for the corrected distribu- 
tions. The mass dependence of the acceptance depends 
on the MC line shape and so systematic uncertainties of 
0.7MeV/c 2 and 2 MeV, respectively are associated with 
the MC line shape choice. These contributions dominate 
all other sources of systematic uncertainty, and the final 
mass and width values are (3914.6i|;|(s£a£) ± 2.0(syst)) 
MeV/c 2 , and (34tg 2 (stoi) ± 5(syst)) MeV, respectively. 

In summary, in the decays B 0,+ —> J /ipwK® ,+ we find 
a J/tpu} mass enhancement at ~ 3.915 GeV/c 2 , confirm- 
ing the BELLE result [2J, but obtain lower mass, smaller 
width, and reduce the uncertainty on each by a factor 
~ 3. The mass is two standard deviations lower than 
the Z(3930) mass, and three standard deviations lower 
than for the AT(3940); the width agrees with the Z(3930) 
and V(3940) values. The ratio of B° and B+ decay to 
YK, Ry, is measured for the first time and found to 
be ~ 3 standard deviations below the isospin expecta- 
tion, but agrees with that for the X(3872) [14]. The 
ratio for the nonresonant contribution Rnr agrees with 
the isospin expectation. We have obtained first measure- 
ments of the branching fractions for B — > J/iptuK and for 
B — > YK, Y — * J/tpoj, for B° and B + decays separately. 
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